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Abstract. At the free-electron laser FLASH, multiple ionization of neon atoms was 
quantitatively investigated at 93.0 eV and 90.5 eV photon energy. For ion charge states 
up to 6-I-, we compare the respective absolute photoionization yields with results from 
a minimal model and an elaborate description. Both approaches are based on rate 
equations and take into acccout a Gaussian spatial intensity distribution of the laser 
beam. From the comparison we conclude, that photoionization up to a charge of 5-1- can 
be described by the minimal model. For higher charges, the experimental ionization 
yields systematically exceed the elaborate rate based prediction. 

PACS numbers: 42.50 Hz, 32.80.Rm, 32.80.Fb, 41.60.Cr 

With the construction of X-ray lasers during the last years applying the free-electron 
laser (FEL) concept of self-amplified spontaneous emission (SASE) [H El [3], H] , a number 
of experiments have become feasible which elucidate the principles of photon-matter 
interaction in a new parameter regime (jSl [B] and refs. therein). The combination of 
high frequency and high intensity allowed in particular new insight into the photoelectric 
effect and the photoionization process [3 El El UHl [HI 1121 1131 [IH US]- Even simple 
photoionization involves in a strict sense multi-electron dynamics as soon as the 
photoelectron is not ejected dominantly from the highest occupied orbital rendering 
a comprehensive theoretical description difficult. For more complex atomic targets such 
as xenon, multiphoton absorption in the energy range of the giant resonance poses 
in addition the question, which role collective electron excitations play in the high 
frequency - high intensity domain [ini HE] . Due to these difficulties, the role of electron 
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correlation for the understanding of multiphoton ionization in the X-ray regime has 
not been settled to date |16l [T71 [18], [191 [20] • Hence, it is worthwhile to carefully 
assess in which cases simple processes, such as sequential photoionization, describe 
the experiments quantitatively and where deviations ask for explanations with a more 
involved dynamics. In the VUV-regime it was found in an experiment with argon atoms 
that sequential ionization describes the measured yields for different charge states well 
^T\. In the hard X-ray regime, it could be shown that sequences of single ionization in 
inner electron shells describe photon-matter interaction in the hard X-ray range and for 
light elements such as neon (Ne) well [12]. The present work determines in a combined 
experimental-theoretical study, to which extend a sequential description for Ne is also 
valid in the EUV regime. For this purpose, results of ion spectroscopy on Ne atoms 
obtained in a focused beam of the free-electron laser FLASH in Hamburg at 93.0 eV 
and 90.5 eV photon energy are compared with a minimal analytical model for sequential 
ionization we have developed and with recent theoretical results from a more elaborated 
approach, including also non-sequential rates [20]. We find that the sequential minimal 
model describes the first five ionization stages reasonably well. The cross sections for 
higher charged ions, where sequential ionization can be excluded, are systematically 
underestimated by the elaborate rate description. Although the effect is weaker as 
compared to Xe with its giant 4d resonance in the EUV [TU], the results point to the 
importance of electron correlation for photoionization, also at higher photon frequencies 
and high intensity. 

The experiments were carried out by applying ion time-of-flight (TOF) spectroscopy 
as shown in Fig. [T] [101 [IH 122] • The FEL beam was focused by the use of a spherical 
Si/Mo narrow-bandwidth multilayer mirror under normal incidence with a focal length of 
20 cm and a refiectance between 60% and 70% [23], depending on the photon energy and 
confirmed with a relative standard uncertainty of 10%. The minimum focus diameter 
as small as rfrwHivi = (2.6 ± 0.5) was derived from the target depletion effect [22] . 
The absolute FEL pulse energy in the /zJ regime was monitored with a relative standard 
uncertainty of 15% on a shot-to-shot basis by means of calibrated gas-monitor detectors 
[23]. The mirror could be moved along the FEL beam in order to shift the focus in back- 
reflection geometry into and out of the interaction volume of our ion TOF spectrometer 
to vary the effective FEL beam cross section [101 H^. As a result, the peak intensity 
(irradiance) of the FEL pulses with pulse durations of AtpwHM = (15 ± 5) fs [16] could 
be varied from 10^^ to 4x10^^ W/cm^ and measured with a standard uncertainty of 
42%. Ne flUed the experimental vacuum chamber homogeneously at the low pressure 
of about 10~^ Pa to avoid any interaction between neighboring atoms and ions. The 
pressure was measured by a calibrated spinning rotor gauge and the temperature by 
a calibrated PtlOO resistance thermometer. The homogeneous electric extraction fleld 
was sufficiently high to collect and register all ions generated within the interaction 
volume of the TOF spectrometer by an open electron multiplier operated in the analog 
mode, i.e., by measuring the charge accumulated on the multiplier anode. The entrance 
aperture of the spectrometer had an extension of 350 /im in the beam direction and 1 
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mm perpendicular to the FEL beam. In front of the experimental chamber, a horizontal 
beam stop of a 1.5 mm in height was introduced which enabled us to diaphragm ions 
produced by the incident unfocused beam (cf. Fig. [T]). TOF spectra were averaged over 
typically 500 consecutive FEL shots with pulse-to-pulse intensity fluctuations varying 
from 35% to 45%. 

Fig. [2] shows typical ion TOF spectra of Ne taken at the highest and lowest intensity 
level, corresponding to the regime of multi- and single-photon excitation, respectively. 
Charged states up to 7+ were observed. Up to 6+, i.e., up to the complete removal 
of the 2p electron shell, the corresponding ionization yields were deduced from the ion 
signals by normalization to the absolute ion detection efficiency and to the number 
of target atoms within the interaction volume. The absolute detection efficiencies for 
singly charged ions were directly measured at low intensity levels where single-photon 
excitation is dominating, using well-known photoionization cross section data of Ne ( [2^ 
and refs. therein). The detection efficiencies for higher charge states were deduced from 
those of singly charged ions assuming that the efficiency is proportional to the ion impact 
velocity [251 126] . The number of target atoms was calculated from the interaction volume 
and the atomic particle density Ua- The latter was determined by Ua = p/kT, where p 
is the target gas pressure, T is the temperature, and k is the Boltzmann constant. The 
interaction volume was calculated taking into account the minimum focus diameter, 
the distance of the interaction volume of the TOF spectrometer from the focus, the 
width of the spectrometer entrance aperture, and the focused beam divergence. The 
latter was obtained using the known focal length of the mirror and measuring the size 
of the unfocused beam incident on the mirror. The relative standard uncertainty for the 
ionization yield is estimated to amount to 36% which arises from the uncertainties for 
target density (2%), interaction volume (30%), and detection efficiency (20%). Fig. |3] 
shows the experimental ionization yields for the Ne charge states from 1+ to 6+ as a 
function of the peak intensity, with filled and open symbols for measurements at 93.0 
eV and 90.5 eV photon energy, respectively. 
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Figure 1. Experimental setup for the investigation of atoms by ion time-of- flight 
(TOF) spectroscopy in a free-electron Laser (FEL) focus of a spherical narrow- 
bandwidth multilayer mirror. 
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Figure 2. Ion timc-of-flight (TOF) spectra of Ne taken at the photon energy of 93.0 
eV and the peak mtensity of (a) 4x10^5 W/cm^ and (b) 2x10^3 w/cm^. In the TOF 
regime below 2.6 /iS, the ion intensities were multiplied by a factor of 30. 



The dependence of the different ionization yields on the peak intensity was also 
calculated by means of a minimal analytical model which is based on the stepwise 
removal of electrons by the sequential absorption of photons, i.e., Ne^ Ne+ . . . 
Ne""*", through a set of rate equations [27] 



iVo(t)= -aoiF(t)iVo(t), 

iVi(t) = aoiF(t)iVo(t) - ai2F(t)iVi(t), 

Mt) = cruF{t)N,{t) - a23F{t)N2{t), 



(la) 
(lb) 
(Ic) 
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Figure 3. Experimental data for Ne ionization yields measured at 93.0 eV (■) and 
90.5 eV (t) photon energy with pulses of 15 fs duration compared with results from 

the present minimal model ( ) and results from the elaborate rate description for 

30 fs pulse length of Ref. [20! ( ). 



with F(t) = Fofit) where Fq is the peak photon flux per unit area and f(t) the 
dimensionless temporal shape of the photon pulses, ajk denotes the single photon 
ionization cross section from charge state j to k. With the new time variable 



r(t) 



dt'fit') 



(2) 



which corresponds for infinite time to an effective pulse length Too = r(t— t-oo), the set 
of Eqs. [T]may be rewritten as: 



dN{I,T) 
dr 



M N(/, 



(3) 



with the time-independent matrix 
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Table 1. Single-photon cross sections ajk in units of 10 ^^cm^, from experiment 
PldH] (bold), using the GIPPER code ^ (italic) and a fitting procedure [30] ■ 

Coi cri2 cr23 cr34 

4.5 ^.07 5.55 5.7i 
4.31 4.55 3.96 



and N(/, r) = {A^'o(/, r), A^i(/, r), A^2(-^, t), • • -j^ where the photon flux Fq is replaced 
by the peak intensity I = Fq hu. The formal solution now reads as: 

N(/,r) =exp(Mr)N(0). (5) 

With the initial condition N(r=0) = N(t-^ - oo) = {1,0,0, ...}^ the vector 
components iV„>o(/,T) may be interpreted as the respective ionization probabilities 
while No{I, r) is the remaining fraction of neutral atoms. 

From this general sequential ionization description we arrive at our minimal model by 
using identical cross sections ajk = cr for all jk, i.e., we need only a single external 
parameter a. For the first ionization steps, the use of identical step-wise cross sections 
can be justified since Uoi, cri2, cr23, and (734 available in the 90 to 95 eV photon energy 
range, do not vary by more than 12% from their mean value of about 4.2 x 10~^^ cm^ 
as can be seen in Tab. [U 

In this minimal model one can solve ^ analytically with the explicit expression 

A" 

Nn{I,Too) = — exp(-A) , A = alToo/hu. (6) 



For comparison with the experimental results, the theoretical data has to be 
convoluted with the spatial intensity distribution of the focused FEL beam within the 
interaction volume. We have chosen 

^(^0' ^' ^) = 1 , /"/ .2 ( ^77^77-^ I ' (^) 

1 + (z/zk) V ^0 [1 + (^/^r) J / 

i.e., a Gaussian distribution in radial direction r. Its standard deviation increases with 



z along the beam from the value Wq = (iFWHM/2\/ln 2/2 = 2.208 /^m at 2; = while 
zn = 0.13 mm denotes the estimated Rayleigh length of the focus for our experiments. 
The results of the minimal model according to Eq. ^ (with a = 4.2 x 10^^*^ cm^. 
Too = A/vr/4 In 2 AtpwHM = 16 fs, hu = (93.0 + 90.5)/2 eV) convoluted with the spatial 
intensity distribution of Eq. ([T]) are shown as a function of the peak intensity / with 
solid lines in Fig. [3] For comparison, also the results of the elaborate rate description 
of the Ne ionization yields are shown by dashed lines pO]. 

As can be seen in Fig. |3| the experimental ion yields from Ne^"*" to Ne^^ are 
well reproduced by our minimal analytical model with the parameter A from Eq. ^ 
globally scaled by a factor 1.20 within the experimental uncertainties for A of about 26%. 
The agreement between experimental data and minimal model within the experimental 
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uncertainties demonstrates both, the sequential character of the first few ionization 
steps, and that the photon intensity distribution of Eq. ^ is appropriate for the 
spatial averaging. The discrepancies between experiment and minimal model at high 
intensity (> 10^^ W/cm^) might be explained by out-of-focus radiation, due to spherical 
aberrations, hitting the neutral atoms which are continuously distributed within the 
vacuum chamber. The predictions of the elaborate rate description from Ref. [20] for 
the Ne^"*", Ne^"*", and Ne^"*" yields are in the same order of magnitude, although a bit 
lower which is most likely due to the fact that the laser pulse length is with 30 fs longer 
than in the experiment (16 fs) but may also be due to the different way of spatial 
averaging. The yield for Ne'^"'" and higher ionization from Ref. |20] , on the other hand, is 
significantly lower than the values of our minimal analytical model and the experiment. 
Note, that the threshold energy for the Ne^"*" — )■ Ne^"*" ionization process amounting 
to about 97 eV [31] and, hence, exceeds our photon energies of 93.0 eV or 90.5 eV, 
respectively. As a consequence, this ionization step starting from Ne^^ requires at least 
two photons, so that rendering its strength much weaker. In our minimal analytical 
model, sequential single-photon ionization is assumed but without explicit reference to 
the ionization thresholds. Interestingly, the yields for Ne'*"^ and Ne^+, which cannot be 
reached by fully sequential photo ionization, are nevertheless well reproduced by the 
minimal model. Starting from Ne^^, the minimal model produces unrealistically large 
yields. 

On the other hand, the data sets in Fig. |3] clearly demonstrate that the predictions of 
the elaborate rate model from [20] are too low compared to the experimental ionization 
yields for the higher charges Ne^^, Ne^"*", and Ne^"*". These discrepancies reveal the 
limitations of rate based approaches, similarly as in the multiphoton ionization of Xe in 
the same photon energy range [IDl [IE] • 

In summary, we have demonstrated that sequential ionization schemes describe 
the generation of lower charge states in multiphoton ionization by soft X-rays well. 
For the higher charges, our experimental ionization yields of Ne are systematically 
underestimated by the elaborate rate description, as already known from previous results 
for Xe. One may conclude that in nonlinear multiphoton-matter interaction electron 
correlation can be significant and should be taken into account to interpret experiments 
at the new X-ray laser facilities. We thank the FLASH team for a very successful FEL 
operation. Support from the Deutsche Forschungsgemeinschaft (DFG RI 804/5-1) is 
gratefully acknowledged. 
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Appendix 

Appendix A. Analytical solution to the coupled rate equations for 
sequential ionization 

In the general case the matrix M in (|5| reads 

/-So ... \ 
+So -Si ... 
M= +si ... 

V ... +Sm-1 / 

where the = o"j^j+iFo are the single-photon absorption rates for taking the 
the charge state i to i + 1. For simplicity we have assumed a maximal charge 
Thus, with the neutral fraction Nq, the matrix M has the size (m+1) x (m+ 
some algebra, one can find a transformation to the diagonal form 

M = T-^MT 

f-so ... 0\ 
-si ... 




\ 








-Sm-1 



with the transformation matrix 
/ 1 



pm- 
^0 






1 

Pl/Q\ 



/ 









m-2 IQ 



P. 

-PT-llQ 



2 

m-2 
m-2 
m-2 









(A.l) 

ion from 
state m. 
1). After 

(A.2a) 
(A.2b) 



0\ 








1 



A.3) 



whereby the first two of the following products have been used 



pi 



n 



Si 



Qk = Yl{s,+i-Sk), 



(A.4a) 



(A.4b) 



The latter two will be used below. 



Note that 



Qi 



I ioi k > I. The 



denominators in the matrix elements of T require that 7^ Si 7^ S2 7^ • • • 7^ Sm-i. 
time evolution of the transformed system N(r) = T^^N(r) is simply given by 



The 



N(r) = exp rM N(0) 



(A.5) 
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In general the initial condition is given by N(0) = T ^N(O). For the particular form 
N(0) = {1, 0,0,.. .y, it simplifies to 

iV„<™(0) = P^-'/R^ (A.6a) 

iV^(O) = 1. (A.6b) 

Finally, we can transform back to the original system by N(r) = TN(r) and obtain 



n-l pn-1 

Nn^„{T) = iV„(T) + ^ -A^NAt) 



exp (— Sj-r) 

Qn 



(A. 7a) 



j=0 

m—l 

NUr) =l-5^iV,(r) 

3=0 

whereby we have used Pq""^ = Po'^Pr^ and jS^ = RiQ] in Eq. (|A.7a|) and J^n ^n{r 



(A.7b) 



1 in Eq. (A. 7b). The two Eqs. (A. 7) give the time-dependent probabilities Nn for the 



various charge states n. 

For the minimal model discussed in the text we need the expressions for degenerate 
absorption rates sq = si = S2 = ■ ■ ■ = s. These can be obtained by subsequent limits 



So — )■ s, Si — )■ s, etc. of Eq. (A. 7a). One obtains the compact form 



ST 



exp (— sr) 



(A.8) 



and Nfn{T) as in Eq. (A. 7b) 
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